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Abstract 

A hybrid experiment consisting of emulsion chambers, burst detectors and the Ti- 
bet II air-shower array was carried out at Yangbajing (4,300 m a.s.l., 606 g/cvo?) 
in Tibet to obtain the energy spectra of primary protons and heliums. From three- 
year operation, these energy spectra are deduced between 10^^ and 10^^ eV by 
triggering the air showers associated with a high energy core and using a neural 
network method in the primary mass separation. The proton spectrum can be ex- 
pressed by a single power-law function with a differential index of— 3.01ib0.11 and 
-3.05 ± 0.12 based on the QGSJET-FHD and SIBYLL-hHD models, respectively, 
which are steeper than that extrapolated from the direct observations of — 2.74ib0.01 
in the energy range below 10^^ eV. The absolute fluxes of protons and heliums are 
derived within 30% systematic errors depending on the hadronic interaction models 
used in Monte Carlo simulation. The result of our experiment suggests that the 
main component responsible for the change of the power index of the all-particle 
spectrum around 3 x 10^^ eV, so-called "knee" , is composed of nuclei heavier than 
helium. This is the first measurement of the differential energy spectra of primary 
protons and heliums by selecting them event by event at the knee energy region. 

Key words: cosmic rays, 7-family, neural network, proton, knee energy region 
PACS: 98.70.Sa, 95.85.Ry, 96.40.De, 96.40.Pq 



* Corresponding author. 

Email address: huaiig@icrr.u-tokyo.ac.jp (J.Huang). 



2 



1 Introduction 



The energy spectrum of cosmic rays is described by a power-law function in 
a wide energy range from about 10^° eV to 10^° eV, however, it shows shght 
changes of the power-law index at several points. These break points of the 
power-law spectrum are assumed to be related to the origin, acceleration mech- 
anism and propagation mechanism of cosmic rays in the Galaxy. One of the 
break points, in which the present paper is concerned, is traditionally referred 
to as the "knee" located around 3 x 10^^ eV. Many experiments have reported 
[1] that the power-law indices below and above the knee approximately take 
the values —2.7 and —3.1, respectively. Although existence of the knee has 
been well established experimentally, there are still controversial arguments 
on its origin. One of them is a possibility that acceleration mechanism could 
be less effective above the knee. Along with this line, there is a general consen- 
sus that stochastic shock acceleration at supernova blast waves could explain 
the cosmic-ray spectrum up to about Z x 10^^ eV [2], or perhaps even higher 
to Z X 10^^ eV [3], where Z denotes the atomic number, despite lack of direct 
evidence. Another argument attributes the knee structure to the leakage of 
the cosmic rays from the galaxy [4] . It is noted that both scenarios mentioned 
above give a rigidity-dependent cutoff for each chemical component leading 
to a heavy-enriched composition of primary cosmic rays at the knee. On the 
other hand, there is another approach in which cosmic rays around and be- 
yond the knee arc assumed to be of extra-galactic origin such as the active 
galactic nuclei [5] or gamma-ray bursts [6]. In this case, the primary chemi- 
cal composition is expected to become proton-enriched. There have been some 
calculations of the primary cosmic-ray energy spectrum based on various mod- 
els on the origin of the knee [7], but all of them are still under debate due to 
lack of detailed knowledge about the chemical composition around the knee. 

Among primary cosmic rays, protons are the key component for understand- 
ing the origin of the knee. Direct measurements of primary cosmic rays on 
board balloons or satellites are the best ways, however, the energy region cov- 
ered by them are limited up to 10^^ eV. The chemical composition of primary 
cosmic rays around the knee, therefore, has been studied with ground-based 
air-shower experiments and/or air Cherenkov telescopes. Since the sensitivity 
to the mass separation among cosmic-ray nuclei with ground-based experi- 
ments is limited, only gross features such as average mass number have been 
discussed. A lot of reports have so far been made on the energy spectrum as 
well as the chemical composition of primary cosmic rays, however, there are 
still serious disagreements among them especially on the chemical composition 
[8]. 

It is possible, however, to improve the sensitivity of an air shower experiment 
to the primary cosmic-ray mass separation by adding a function to observe 
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the energy-flow characteristics of air-shower cores at a high-mountain altitude. 
The Tibet hybrid experiment [9] was designed to detect 7-families in an air 
shower core by large-area emulsion chambers in coincidence with an accom- 
panied air shower, where a bundle of energetic 7-rays and electrons detected 
by the emulsion chambers are called a 7-family, which is caused by a young 
air shower. Among primary cosmic rays, protons and heliums can penetrate 
deep into the atmosphere and produce a young air shower accompanied by 
a 7-family event most efficiently due to their longest interaction mean free 
paths in the air. Therefore, tagging an air shower with a 7-family event en- 
riches the proton and helium component naturally. Another merit in doing a 
hybrid experiment in Tibet is that the atmospheric depth of the experimental 
site (4300 m a.s.l., 606 g/cm^) is close to the maximum development of the 
air showers with energies around the knee almost irrespective of the masses of 
primary cosmic rays. We can determine the primary cosmic-ray energy much 
less dependently upon the chemical composition [10] than those experiments 
at the sea level. Thus, the Tibet hybrid experiment enables us to measure the 
primary proton and helium differential energy spectra by selecting them event 
by event. 

In this paper, we briefly report on the study of the energy spectra of the 
proton and helium components in cosmic rays around the knee energy region 
obtained with the Tibet hybrid experiment. 



2 Experiment 

The Tibet hybrid experiment, consisting of emulsion chambers (ECs), burst 
detectors (BDs) and the Tibet-11 air-shower array (AS), composed of 221 
scintillation counters each placed on a 15 m square grid with an enclosed area 
of 36,900 m^ was operated at Yangbajing in Tibet during the period from 
November 1996 through August 1999 [9] and a total live time of 699.2 days. 

The AS is used to measure the shower size and the arrival direction of each air 
shower. Any four-fold coincidence in the detectors is used as the trigger con- 
dition for air-shower events. The air shower direction can be estimated with 
an error smaller than 1°. The primary energy of each event is determined by 
the shower size (Nc). The energy resolution is estimated to be 17% at ener- 
gies around 10^^ eV by our simulation, almost independent of the interaction 
models used. 

The ECs and the BDs are constructed near the center of the AS [11], and 
are used to detect high-energy air shower cores accompanied by air showers 
induced by primary cosmic rays with energies above ~ 10^"^ eV. The total area 
of ECs is 80 m^. The basic structure of each EC is a multilayered sandwich 
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of lead plates and X-ray films (FUJI X-RAY FILM TYPE 200) of 40 cm x 
50 cm in area, where X-ray films are put every 1.0 cm of lead in the chamber. 
The X-ray films in ECs are replaced by new ones every year to reduce the 
background. The ECs are used to measure the energy, the position and the 
arrival direction of each 7-family shower with energies above 1 TeV. The BD 
with the same area are placed just below 4 ECs, namely, 4 ECs are set above 
one unit of the BD. Thus, 400 blocks of ECs and 100 BDs in total are used 
in this experiment. A burst event is triggered when any two-fold coincidence 
of signals from four photodiodes of each BD appears. When the BDs trigger 
an event, its accompanying air shower is simultaneously recorded. The BDs 
are used to measure the burst size iVb and the position of each air shower 
core. The arrival direction of the 7-family event is determined by the spatial 
reconstruction of the cascade showers in ECs, whose details are described in 
[12]. The matching between an AS and a BD event is made by their arrival 
time stamps, and the matching between a 7-family event in EC and the BD 
event is made by their positional correlation, and the matching between the 
7-family event in EC and the AS event is made by their directional correlation. 

In the following analysis, we present our results based on the ECs and AS 
data, as those obtained from the BDs data were published in the previous 
paper [11]. 



3 Simulation 

We have carried out a detailed Monte Carlo (MC) simulation of air showers 
and 7-familics using the simulation code CORSIKA (version 6.030) including 
QGSJETOl and SIBYLL2.1 hadronic interaction models [13]. From a point of 
view to check the dependence of the obtained results on the assumed primary 
cosmic-ray composition in MC, two primary composition models are examined 
as the input energy spectra, namely a heavy dominant (HD) and a proton 
dominant (PD) ones [9]. The energy spectrum of each component in the HD 
model has a rigidity-dependent break point of the power index with proton's 
knee around 1.5 x 10^^ eV leading to the dominance of the heavy component 
at the knee energy region. On the other hand, it is assumed in the PD model 
that light components are dominant up to the knee, in which every component 
has the same break point of the power index at the knee energy. In both 
models, the fraction and the power index of each component are determined 
by fitting to the fiuxes of the elements obtained by direct observations below 
10^^ eV, and fitting the sum of the each element at higher energies to the 
all particle flux obtained by air shower experiments. Therefore, the difference 
between two models exists in the fraction of the elements above 10^"^ eV. The 
fractional contents of the assumed primary cosmic-ray flux models are listed 
in Table 1, together with those for making air showers accompanied by 7- 
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Table 1 

Fractions of the proton(P), helium(He), medium(M) and iron(Fe) components in 
the assumed primary cosmic-ray spectrum of the HD and PD models(upper table), 
together with those for making air showers accompanied by 7-families (lower table) 
(see the Section 4). 
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families, where M denotes the sum of medium heavy elements between helium 
and iron. One can see from Table 1 that 100% of the 7-family events are 
induced by protons and heliums below 10i5 eV. At the higher energy range, 
however, the contribution of the other nuclei heavier than helium increases 
with the energy which amounts to around 15% and 4% in case of HD model 
and PD model, respectively. The method of the separation of primary mass 
groups using a neural network method is described in the next section. 

The numbers of generated 7-family events satisfying the criteria described 
in the next section are 5252, 5926, 8588 and 7376 for the QGSJET-t-HD, 
QGSJET+PD, SIBYLL+HD and SIBYLL+PD models, respectively, which 
are more than 30 times greater than the experimental statistics. 



4 Analysis 

Shower spots registered by high-energy 7-rays or electrons in the X-ray films 
were automatically analyzed by using image scanners [12]. The 7-family events 
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are selected by imposing the conditions of Ej^^ = 4 TeV, A^-^ > 4, ^ E^'^ > 20 

i 

TeV and < R >> 0.2 cm, where E'^*'* is the minimum energy for a cascade 
shower, the number of cascade showers in a 7- family, J2 the sum energy 
of cascade showers in a 7-family and < R > {< R >— y^r^'/A^^) the mean 

i 

lateral spread in a 7-family. In this experiment, we observed a total of 177 7- 
family events, each of which is accompanied by an air shower with Ne > 2x 10^ 
corresponding approximately to 5 x 10^^ eV for a proton. 

The separation of the primary mass is made with use of a feed-forward artifi- 
cial neural network (ANN [14]) whose apphcability to our experiment was well 
confirmed by the Monte Carlo simulation [15,11]. The first task of the ANN 
is to separate protons from everything else by training the network with a 
proton flag. As shown in Table 1, heliums are the main component of the con- 
tamination at lower energies while medium heavy component also contributes 
in higher energies depending on the primary composition. Then the second 
task of ANN is to separate a group of proton+helium from others above 10^^ 
eV by training the network with a light element flag, i.e., a flag for proton or 
helium. Thus we can obtain helium spectrum by subtracting the number of 
proton events obtained in the first task from the proton-|-helium dataset. 

We examine four cases, QGSJET + HD, QGSJET + PD, SIBYLL + HD and 
SIBYLL + PD models. The following six parameters are input to the ANN 
with 30 hidden nodes and 1 output unit, which are abbreviated to a 6 : 30 : 1 
network ^ : (1) E-E-y in ECs, (2) A^^ in ECs, (3) < R > in ECs, (4) the mean 
energy-flow spread < ER >— '^E^^r-y'^/N^ in ECs, where Ej and r-y* are 

i 

the energy of each cascade shower in the 7-family and its distance from the 
energy-weighted family-center, respectively, (5) the air shower size A^e in AS, 
and (6) sec^, where 9 is the zenith angle of the arrival direction. 

One MC data set was generated for each of the four models and divided 
into two subsets; the one is used for training the ANN, and the other for 
estimating the ability of the ANN to classify the nuclear species. Then, the 
training MC data subset is fed to the ANN in a number of training cycles 
of 2000. ^ To train the ANN in separating protons from others, the input 
patterns for protons and other nuclei are set to and 1, respectively. After 
the training, the another MC data subset is used to estimate the purity and 
the selection efficiency of protons. Then, the ANN output pattern value (T) 
is a real number from to 1. The T distributions in the QGSJET+HD model 
is presented in Fig. 1, together with the experimental data. One can see that 
the experimental data is in a good agreement with the MC prediction, and 

^ We compared the results from various number of hidden nodes and found the 
results are not changed when we use the number of hidden nodes greater than 20. 
^ The training curve becomes stable at number of training cycles 500. 
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Fig. 1. ANN output pattern value (T) distributions compared with MC 
(QGSJET+HD model). Used symbols are ; experimental data (the closed squares), 
MC:all (the solid-red histogram), MC:protons (the dash-dotted green line), 
MC:other nuclei (the dotted-black line). 



that the proton-induced events are clearly separated from other nuclei. ^ We 
define a critical value of T to separate protons from others requiring the high 
purity, which reduces the effect of the contamination, and the high selection 
efficiency of the proton events, which reduces the statistical error. Since these 
two factors are of the competing nature, the purity and the selection efficiency 
are calculated as a function of the critical value Tc and its value is set as 
0.4 where average purity and selection efficiency over whole energy range are 
greater than approximately 85% and 70% for all models, respectively. It was 
confirmed that the final result does not depend on the different choice of the 
Tc around 0.4. The high value of the purity as listed in Table 2 and their 
mutual deviation among different models in the same energy interval being 
within 4% assures us the quality of the proton selection. 

Finally, we obtained 111, 111, 112, 112 proton-like events out of 177 observed 
events after the 7-family selection based on the QGSJET+HD, QGSJET+PD, 
SIBYLL+HD and SIBYLL+PD models, respectively, among which 110 events 
are identical and one event belongs only to QGSJET analysis and two events 
belong only to SIBYLL analysis. 



3 We also examined the QGSJET+PD, SIBYLL+HD and SIBYLL+PD models, 
obtaining similar results to the ones as shown in Fig. 1, for the selected events 
are already mostly proton and helium origin in all models as shown in Table 1. 
Of course there are some differences between HD model and PD model, reflecting 
the degree of the contamination, however, it is not large enough to rule out one of 
them because of experimental statistical errors. This figure demonstrates how the 
proton-induced events are separated from those induced by other nuclei. 
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Table 2 

The purity of the selected events by T < 0.4. 



Model 


Energy (eV) 


Purity(%) 






HD 


PD 


QGSJET 




96.7±0.7 


97.4±0.4 




lOis.iQie 


83.1±1.6 


86.7±0.8 


SIBYLL 




96.2±0.5 


97.3±0.3 




lOis.iQie 


82.8±1.2 


86.1±0.7 



5 Results and Discussions 



In Fig. 2, we show the measured primary cosmic-ray proton energy spectra 
assuming the two interaction models (QGSJET and SIBYLL) and two pri- 
mary composition models (HD and PD), together with the results from other 
experiments. As seen in Fig. 2, the present results assuming the HD and PD 
models in the simulation are in a good agreement with each other within the 
statistical errors. The measured proton energy spectra can be expressed by a 
single power-law function of a differential form J(£')(m~^-s^^-sr~^-GeV~^) = 
A X 10-13 X (_^)-B^ where {A,B) is (4.56±0.46, 3.01±0.11), (4.14±0.44, 
3.08 ±0.11), (3.21 ± 0.34, 3.05 ± 0.12) and (3.24 ± 0.34, 3.08 ± 0.12) based on 
the QGSJET+HD, QGSJET+PD, SIBYLL+HD and SIBYLL+PD models, 
respectively, where the errors quoted are the statistical ones. The error in the 
spectral index is statistics dominant, while that in the absolute flux value is 
model-dependence dominant. For the absolute flux value, the QGSJET model 
gives approximately 30% higher flux than the SIBYLL model. This can be 
mainly attributed to the difference of Feynman x^-distribution of charged 
mesons between QGSJET and SIBYLL model in the very forward region at 
a collision [13]. The Feynman a^i^-distribution in the SIBYLL model is harder 
than that in the QGSJET model in the xp > 0.2 region, so that the gener- 
ation efficiency of 7-families by the former model becomes higher than the 
latter, resulting in a lower proton flux in the case of the SIBYLL model. As 
compared in Fig. 2, the present results are consistent with those obtained by 
the burst detectors in this experiment within 25% [11]. This imphes that the 
systematic energy-scale uncertainty in our experiment is estimated to be 10% 
level. A solid straight line with the power index —2.74 drawn in Fig. 2 is the 
best fltted line for the data points in the energy region below lO^"' eV observed 
by recent direct measurements [22], which is harder than the indices of our 
proton spectra. 

Thanks to its light mass, the helium component can also trigger our hybrid 
experiment although the efficiency at 10^^ eV is about 4 times lower than the 
case of protons. The ANN method is again apphed to obtain the helium spec- 
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Fig. 2. Energy spectra of primary cosmic-ray protons obtained by the present ex- 
periment (a) and they are compared with other experiments (b): Tibet-B.D.[9], 
KASCADE[16], JACEE[17] and RUNJ0B[18]. The all-particle spectra are from the 
experiments : PROTON satellite [19], Tibet-III[20] and AKEN0[21]. For the solid 
line with the power index —2.74, see the text. 
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Fig. 3. Energy spectra of primary cosmic-ray helium nuclei obtained by the present 
experiment (a) and they are compared with other experiments (b) 

trum over the energy 10^^ eV. Because of the training algorithm of ANN, it is 
not possible to train the network to separate heliums from others directly, for 
the helium mass is between protons and other heavy nuclei and the character- 
istics of the helium event is smeared out by the fluctuation tail from the both 
sides. Therefore we train the network to separate light component (proton or 
helium) from other nuclei, by assigning to light component and 1 to other 
nuclei. The critical value to select light component is set as 0.2 where the 
selection efficiency reaches to 70% and the purity is 93% for all models. Then, 
the helium spectra can be obtained by subtracting the number of protons, 
which are previously obtained by proton-training, from the number of pro- 
ton-|-helium events. Above mentioned procedure was applied on each energy 
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Fig. 4. Fraction of the primary cosmic-rays heavier than hehum nuclei obtained by 
assuming the QGSJET and SIBYLL interaction models. Our results are compared 
with those by the KASCADE experiment [16]. 



bin to obtain the energy spectra of heliums and the result is shown in Fig. 3, 
where the same dependence of the absolute intensity on the interaction models 
is seen as in the case of proton spectra. 



We can also estimate the fraction of the nuclei heavier than helium in cosmic 
rays around the knee using the proton+helium spectra and the all-particle 
energy spectrum obtained by the Tibet air shower array [20]. Shown in Fig. 4 
is the fraction of primary cosmic rays heavier than helium nuclei assuming 
the QGSJET model and the SIBYLL model which are compared with those 
obtained recently by the KASCADE experiment [16]. Our results using 4 kinds 
of simulation models commonly indicate the average mass of primary cosmic 
rays is going up around the knee, towards the direction of heavy dominance. 
On the other hand, the KASCADE results which measures both air shower 
size (Ne) and muon size (A^^) to deduce the energy spectrum of separate 
mass groups from the all-particle energy spectrum, strongly depend on the 
interaction models. The muon size contained in the air shower depends on the 
number of charged pions produced in the central and backward region (in the 
center of mass system) in the collisions of primary cosmic rays on air nuclei, 
which has a sizeable uncertainties experimentally as well as theoretically and 
is largely dependent on the interaction models. From this point of view, the 
size of low-energy muons may not be a suitable parameter for separating 
the air showers into different primary mass groups. 
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6 Summciry 



A hybrid experiment of emulsion chamber and air-shower array was success- 
fully done at Yangbajing in Tibet to study the primary cosmic rays around 
the knee energy region. Using the events observed simultaneously in the emul- 
sion chamber and the air-shower array, and applying a neural network analysis 
to this data set, we obtained the energy spectrum of primary protons in the 
energy range from 4 xlO^^ eV to 10^^ eV. The spectrum observed can be 
represented by the power-law fit and the power indexes are estimated to be 
—3.01 ±0.11 and — 3.05±0.12 for the spectra obtained using the ANN trained 
by the QGS JET+HD and SIBYLL+HD events, respectively, which are steeper 
than that extrapolated from the direct observations of —2.74 ±0.01 in the en- 
ergy range below 10^^ eV. The absolute fluxes of protons was derived within 
30% systematic errors depending on the hadronic interaction models adopted 
in the Monte Carlo simulation. We also estimated the primary helium spec- 
trum at energies above 10^^ eV, which has almost same spectral slope with 
the proton spectrum. 

We further obtained the result that the fraction of the nuclei heavier than 
helium in the primary cosmic rays around the knee region, which was estimated 
using the proton±helium spectrum and the all-particle spectrum observed 
with the Tibet experiment, increases with increasing primary energy. This 
strongly suggests that the main component responsible for making the knee 
structure in the all-particle energy spectrum is the nuclei heavier than the 
helium component. 

This is the first measurement of the differential energy spectra of primary 
protons and heliums by selecting them event by event. In the very near future, 
we will start a new high-statistics hybrid experiment in Tibet to clarify the 
main component of cosmic rays at the knee. The new experiment is able 
to observe the air shower cores induced by heavy components around and 
beyond the knee, where direct measurements are inaccessible because of their 
extremely low fluxes [23]. 
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